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We present a method for determining a purely quark Lagrangian by mocking up the QCD partition 
function for large gauge couplings g. The resulting effective theory displays all the symmetries of 
low energy QCD and can be potentially used to explore hadron properties. 
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I. INTRODUCTION 


Quantum chromodynamics is an SU (3) gauge theory which contains quarks and gluons that interact with each 
other. Whereas at high energies the coupling constant is small and one can adequately describe the reality by 
expanding perturbatively in the coupling constant it is still not completely clear how one can depict the low-energy 
regime where the coupling constant is large and quarks form bound states of mesons and baryons. Various methods 
and Lagrangians have been proposed to describe the low energy dynamics: While some of them like the Nambu Jona 
Lasinio model [Il:i are still based on fermions as fundamental degrees of freedom, most of them consider the mesons 
and baryons as the starting point. In the latter case one then constructs effective low-energy models which are widely 
based on the symmetries of the QCD Lagrangian like the SU{3)l x SU{3)n chiral symmetry, U{l)v or U{1)a (as 
well as assumptions about QCD vacuum and potential) embedded in the quark flavor sector of QCD. Models of low 
energy QCD with notable results regardiM hadrons properties and interactions include chiral perturbation theory 
3,(3], linear and nonlinear sigma models U-[il| or other symmetry induced effective Lagrangians [2^- 321. Thus one 


can conclude that chiral symmetry whether spontaneously or explicitly broken is the major ingredient for building a 
low energy QCD effective theory. 

In the present work however we shall consider a completely different point of view: we will not rely on symmetries 
but instead we shall build a purely dynamical Lagrangian and obtain the chiral symmetry as a derivative of the 
method. Our main tool is the QCD partition function from which one can derive all the properties of the particles 
and interactions. We start by considering an alternative description of the partition function after integrating out 
the quark degrees of freedom. As a consequence the original quarks will be replaced by copies (that describe the 
constituent quarks) with the same masses and quantum numbers but with different Lagrangian and thus interactions. 
In the end we shall obtain a Lagrangian where the gluon degrees of freedom have been eliminated and that contains 
only the quarks and their subsequent interactions. Note that our method implies only manipulation and calculations 
of the partition function and does not involve any loop computation. 

We consider the effective fermion Lagrangian obtained in this paper as a purely theoretical one and as a first step 
towards a more comprehensive approach. The connection with the phenomenology of the bound states of mesons and 
baryons will be performed and finalized in a future work. 


II. THE QCD LAGRANGIAN AND PARTITION FUNGTION 

We start with the gauge fixed QCD Lagrangian with N colors and Nf flavors in the fundamental representation: 

^ + E - rufYf, (1) 

/ 


where, 


F^,=dY:-d.Al+gr’>^AlAl 


( 2 ) 
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and, 

D^ = d^- igApr (3) 

Here as usual are the generators of the group SU{N) in the fundamental representation. 

We shall ignore the ghosts in what follows as they do not contribute essentially to our arguments. First we separate 
the Lagrangian in two pieces: 

= -ruf + 

/ 

^2 = (4) 

The QCD partition function has the form: 

Zq = J dA‘lj^{x)d^{x)dc'^{x)d^fixd'i/fi{x)exp[i j d'^xC], (5) 

where the usual rules of the path integration apply. Here / is the flavor index whereas I is the color one. Then one 

can integrate over the fermion variables to obtain: 

Zo = y (x) det [i'y^df, -ruf + 57 '"^“^“] exp [J d'^xC 2 ]- ( 6 ) 


III. A SIMPLIFIED APPROACH 


We shall discuss in particular the real life case of iV = 3, Nf = 3 corresponding to QCD with three light flavors. 
We introduce the fermionic current that couples with the gluon field for each fermion species: 


j;f = ( 7 ) 

For one generation of fermions there are 8x3 degrees of freedom for each fermion where 8 represents space time 
degrees of freedom of an off-shell fermion and 3 the number of colors. Overall we have 24iV/ degrees of freedom for 
all fermions. 

We consider a single flavor and next extend our arguments to Nf flavors. We start by making a change of variables 
from the elementary fermion degrees of freedom to the composite current J“. We first mention that in the case of 
anticommuting variables the Jacobian appear with an inverse power as compared to the case of commuting ones. Note 
that for this change of variables to make sense we need 24 degrees of freedom for J“ instead of the 32 that one might 
obtain from a simple counting. This means that we need something similar to a gauge conditon which we choose to 
be J“ = uj where uj is an arbitrary function as in the more standard case of a gauge field. Since there are exactly 8 
constraints one obtains the desired matching of 24 degrees of freedom for both the fermion and vector boson variables. 
Then the change of variable is: 


J d^d'i) J dr^s (S'" j“ - uj) 


dJ?. 


d'^r 


( 8 ) 


The Jacobian in Eq. m is a determinant with dimension 24. Since any arbitrary derivative 75 ^ = jminti) 

(here the first indices are space time whereas the second ones are color) contains a fermion variable then the deter¬ 
minant will contain products of 24 fermion variables. Taking into account that there are exactly 24 distinct fermion 
variables and the anticommuting nature of these we conclude that the actual determinant will be given exactly by the 
product of the 24 distinct variables (since those term that contain a repeated variable will be zero) times an irrelevant 
constant factor. Thus: 

= const n ^ im • (9) 

i,m 

Then one can also write quite safely: 

^ detrrj^] ( 10 ) 
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Here the determinant is taken in the space so it has the dimension 12 which leads exactly to a product of 24 
distinct fermion variables. Note that at each point we take into account the composite nature of J“. 

We need to extend our arguments to Nf = 3 flavors. By applying the previous approach (or simply consider a 
product over the umber of flavors we obtain: 


Y[d^fd^f 

/ 




n dr,f deti^er^f]s - u) 

f 


(11) 


It is more convenient however to write: 


ndet[7^t“J“/] 


const X det[t“7'^^ 


/ 


( 12 ) 


The above relation is due to the fact that the above determinant contains a product of 72 fermion variables. But this 
is the total number of fermion variables so that each term that contains a variable twice is zero. Thus the determinant 
will be a constant times the product of 72 distinct fermion variables. But this is also what one would obtain from Eq. 
(HU so Eq. (HU is correct. 

This being settled we need to write the full partition function in Eq. m in terms of the new variables J“j. As 
we mentioned previously the change of variables makes sense only with the additional constraint on the field J“. We 
shall consider a particular case of it with w(x) = 0. On the other side this constraint is equivalent to constraining the 
free equation of motion of the fermion field. Thus it can be fulfilled only if, 

— m'l' = ^^-4'(x). (13) 

2 / \ 

We use the function (where M is an arbitrary scale) instead of simply a(x) because we do not want to further 

constrain the fermion fields 41. Then the kinetic term for the fermion field will be replaced in the lagrangian by: 

- m)4- = (14) 

Furthermore since the function a{x) is arbitrary one can use the equation of motion to eliminate the kinetic term 
altogether from the Lagrangian. 

The partition function in Eq. (jS]) will thus become: 




= / 

/ 


det[r7' 


! 


exp 



(15) 


The same partition function is obtained if one introduces three flavors of fermion copies %/, Xf of fh® original 
quarks such that: 


Zi= i dA“P[dJ^/, P[dx/dx/exp 




exp 


* [ d^x[Y,r^f,gA'^>^ + C,] 

J f, 


(16) 


We make the change of variable if®! = X)/ K ^2 = to determine: 


[ dAidK;,Kl,dK^,,Y[dxfdxfexp 


X exp 

* f d^x[gK;,A^f^ + C 2 ] 

' f 

. f 


J 




( 17 ) 
















4 


where the integral over K'^-^ is a delta function. Here one can drop the unwanted integrals over since 

they do not contribute in any process and apply the delta function to obtain that the effective Lagrangian is just: 

= C, = -^ E X/7"i“x/j . (18) 

Although this procedure seems oversimplifying it gives a correct glimpse of what kind of Lagrangian we should 
expect in terms of the current J“ if the fermion kinetic term is neglected. The next section will contain a more 
general and comprehensive approach. 


IV. A COMPREHENSIVE APPROACH 

We start with Eq. (jS]) which we rewrite here for completeness: 


Zo = J det [ 17 ^ 9 ^ — m/+ exp i j d'^xC 2 


(19) 


/ 

We shall now try to reproduce the above partition function using a different set of variables. 
We first need an identity that we shall prove in what follows. Consider the integral: 


Z:r) — 


2 ^ 2 _^ 

dx^dxidyjdyjdJkdSm e:x.-p[i[-^xA’^xj'^ + —^yB^yS^y + ^ E 


( 20 ) 


Here yi and their conjugates are each a set of n Grassmann variables which we shall assume have the mass 
dimension and Jk, Sk bare a set of regular commuting ones each with mass dimension m. The index k goes also 
from 1 to m and the matrices are each a set of m (n x n) matrices where have mass dimension m and 

have mass dimension mP. We will solve first the integral over the fermion fields to get: 

Z,= f dj'=d.5'=det[z^A'^j'=]det[*^i?-5™]exp[*^^j'^5"]« 

k 

J dJ'^'dS''"det[^A''H™j'=5'™]exp[i^^j''5''=] « 

k 

fc 

J dS-dzdz. + « 

k 

/ dz^dzi ■ 


, l + ^zA>^B>^z 


f dz^dzi exp[- ^ ln[l + -^zA'^B’^z]] Ki 
f dzidzi exp[— ^ w det[A^H^] 


( 21 ) 


Here we dropped all unimportant constant factors and in the last line we consider only terms in order and dropped 
the higher order ones in the fermion fields. This is the only approximation involved in Eq. (EH). 

Now consider again Z^ and this time we integrate first over and Sk- 

Za;= f dxidxidyjdyjdJkdSm exp[i[^xA'"xj'‘ + ■^baryB'^yS'^ + ^ ~ 

= f dxdxidyjdyjdJk S{JP + -^yB^y) exp[i^xA'=xj'^] ss 

J dxidxidyjdyj e:>cp[—ij^xA'^xyB^y] « det[A'"H'"]. 


( 22 ) 
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The result in the last line of Eq. (12^ will be of most importance in what follows. 

There is an important extension to the results in Eqs. (I^Tl) and which we will use here but state without proof 
because the proof is just a simple generalization of the arguments above: 

Z^= f dxidxidyjdyjdJkdSm - -^yB^yS^ + X! ^ 

= J dxidxidyjdyj exp[i^j[xA^xyB^y]] Ki det[B"^A^]. (23) 

Here C is an n x n matrix but one would obtain the same result if C were a polynomial with terms of the type {xRx)^ 
with R any n x n matrix and k an arbitrary integer. The reason stems from matching the degrees of freedom of 
Grassmann variables x and y. 

We thus start from: 

Zx = J dxidxidyjdyjdJkdSmexTp[i J d'^x[^xA'^xJ^ — yB^yS^ + +T[x, a:)]] = 

k 

= J dxidxidyjdyj exp[i J d^x[^^xA'^xyB'^y + T{x,x)]] det[B'^A’^], (24) 

and consider the set Xi as being fermions 'i’if , where N is the number of colors and is the number of flavors. 
Similarly the set yi is a similar set of fermions Xif- total we consider to have Nf = 6 flavors of fermions 

corresponding to the 6 flavors of quarks of the standard. We further define the matrices A^, B^ as: 

A’" = 

B^ = 




(25) 

(26) 


where n is a constant with mass dimension 1 and mo is an arbitrary scale. Then: 

B^A’^ == - Amee + g'-i^eA^] = 

- 1 - 1 

o Ar -4m -h = 


^ 2N 

- 1 


2N 


2N 


(27) 


where we redefined g' = 9 observe that the operator B^A^ is the operator that appears in 

the standard model between two fermion states. Note that we can include in each and B^ a diagonal flavor matrix 
with the same final result. 

Then the counterpart of the Eq. (|M1) in terms of the above definition will be: 


Zq = J dA“(a;) det — m/+ 57 ^A“] exp i[J d'^xC 2 + J d^xT] 


= / dA“ dT/d^-dx/rfx/exp 


/=i 


d'^xjCo 


exp 


^M3 


3 3 

- mt“7^ + 5'^“)^-/] ^ [x/'7'"i“x/'] + f d'^xT] 
/=! /' = ! d 


(28) 


Here M is an arbitrary scale that reestablishes the correct dimensionality. 
We denote: 


Zi= Y\_dAifdAifdxfdxfexp 
d /=i 


d f=i //=i d 


,(29) 














6 


and work only with it. 

First we will make the change of variable (see section II for details). Note that this implies a 

subsequent gauge condition for which we shall discuss later. The partition function Z\ will become: 


^3 ^ r ^ ^ r 

Zi= d^/d^'/dJ“^det[J“^7''t“]exp i[— / ^ + / d^xT] 

J /=! ■> /=! /=! 


.(30) 


Noting that in the change of variable we can use det[X)/ instead of H/ det[J“j 7 '^t“]^ (see Eq. (fT^ we can 

further write: 


/=1 / 


i[^ [ d^x'^['i>f{it‘^df, - ^ + f d'^xT] 

J f=i f=i J 


(31) 


We can further make a change of variables — '^^2’ ^^3 ~ drop the unwanted integrals 

over and as they would not contribute to any process. This yields: 


dr; / nd^'/dvi/;det[r;7'^t“]3exp 


/=i 


j dl^xT] 

J J 


/ n d'i/fd'i'fd^fd^f 


/=i 


exp 


3 3 

/ d^^^[^fizt‘^d,-mex, + 9 'A‘;,)^f]Y‘^^ + + j d^xT] 

d /=! /=! d 


(32) 


Here we introduced another set and of 3 fermions to account for the determinant in the first line of Eq. 

Here a is an arbitrary dimensionless coupling constant 

As we mentioned previously in order to be able to make a change of variable from the fermions X/ to the currents 
one would need something similar a to gauge condition to cut off the number of degrees of freedom from 32 to 
24. We shall consider this constraint as 9^ J“ = 0. Upon the change of variables to Y;, and C/; this constraint 
will become d^Y^ = 0 for the only variable of interest. We shall introduce this in the partition function Zi as: 


S{d>^Yp = J dS'“exp[i J = J dS'“exp[-i J J“], 


(33) 


where Mi is an arbitrary constant with mass dimension 1. With the addition of the gauge condition the partition 
function in Eq. (15^ will become: 

3 3 

Zi J' dY^ dS /n d^/dvl/^det[Y;Vt“]3exp *[T J d^x^[^fizt‘^d^-mt‘^x^+g'Al)^f]Y‘^^ + J d^xT] 


/=i 


dY;^ / n 


/=i 


exp 


3 3 

* / + 9'Al)^f]Y‘^>^ + a E “ Mid^S<^Y<^^ + T] 

/=i /=i 


(34) 


We shall integrate 1"“^ by observing that it couples only linearly that the Lagrangian is hermitian and thus it leads 
to a product of delta functions. This yields: 


= 




/=i 


S{g'Al E ^f'^f + + MaJ2^n^nf - MiMd^S^) x 


exp[i J d'^xT] 


/=i 


/=i 


(35) 
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We shall regard the delta function in the space of variables . 

Before going further we need to make an important ammendment. When we introduce the change of variable from 
Xf, Xf to we replaced (we consider here for simplicity only one fermion species): 

J dxfdxf- J (36) 

The reason stems form the fact that when we transform from one set of variables one gets always products of 24 
different fermion components. However one can extend the above transformation to include some function in the 
determinant. Consider that we have: 


J dJ^ det[ix^dfj, —m + (37) 

and we apply the inverse transformation to fermion variables %/: 

J dJ^ dei[ix^d^ —m + = 

/ dJ°- 

dxfdxf\^\-^det[trd^-m + krrxn^.t‘^Xf] (38) 

Moreover we also need to consider a function (for one flavor) to be integrated which is of the type det[X“J“]. 
Assume we express this in terms of fermion components. Those variables that repeat themselves will get canceled. In 
the end det[A“J“] = fxiX 2 --XiX 2 where f is an arbitrary function and the product is over all 24 fermion components 
for one flavor (note that the product might contain derivative which we omit for simplicity). Then Eq. (1381) will 
become: 


dj; det[i7'"5;, - m + det[A“ j;] = 

1 


dxdx- 


XlX 2 --XlX 2 - 


[ao + aiXiXj + ■■■ai2XiX2-XiX2-] x [xiX2---XiX2--] = 


J dxdxiao + aiXzXj + ■.■ai2[xiX2---XiX2--]] = j dxdx[ai2XiX2-XiX2-] 


(39) 


But the end result corresponds to that part of det[i7''9^ — m + kt°‘x^J'^f\ that contains only J“ so the addition we 
make to det[t“ 7 ^ J“] is irrelevant. The main point of the above discussion is that even when we integrate over J“ the 
intrinsic fermion nature of the variables should be considered. For most of the purposes we shall consider J“ however 
as a regular commuting variable. However our discussion has a counterpart by discussing purely the J^’s. This being 
settled one can add in Eq. dSSl) a kinetic term for the fermions as in: 


= 


J]^ d^fdd/fd^fd^fdS^ 


/=i 




/=i 


/=i 


exp[i[ J - Mf)^f] + J d'^xT]]. 


(40) 


The expression in Eq. (1401) is the final partition function we will work with. 

It appears that the delta function in Eq. (HDl) is badly defined. In order to fix this we go back to the full Lagrangian 
of interest (including the pure gluon one) which is: 


c = Y,[^f{trd^-Mf)^f]+T+ 

f 

3 3 

+ Ma ^ - MiM5^5“)] + £ 2 (A“, g) (41) 

/ /=i /=i 



We first recall that g' = " 2 N 3 make the change of variables => Note that we can do this at any time 
without affecting in any way the dynamics. This yields (we use for simplicity the same notation): 

f 

— 1_^_ ^ - 1 

E (42) 


/ 


/=i 


/=i 


where £3 does not contain any more any trace of g. We are interested in the regime where g is large so we can consider 
^ a small parameter. Next we take into account that one can solve from the delta function as function of the 
other variables. With this substitution the Lagrangian becomes: 

£ = - Mf)0] +T+ ^f, ^/, 4//, 5“) 


/ 

Assume we consider separately shifts in the variables f and their subsequent conjugates: 

0=^f + 

0 = 0 + 

1 

3'^/ 


(43) 
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^f = 0 + 


(44) 


Then one can still obtain the interaction Lagrangian of order provided that the fields 'k f and satisfy the equation 

of motion for the free Lagrangian Cq. Recalling the definition of T from Eq. (HID we get: 


It = -^4'(V^/^'/ -u^) = 0 

c^'E' run ^ ^ ■' 


(45) 


whose only non trivial solution isX]/^/^'/ = Thus the unwanted denominator in the (5 function in Eq. (1401) is 
conveniently fixed. 

We redehne and ■^MMi = where y and z are two adimensional coefficients. 

From Eqs. (gQl) and we obtain the effective fermion Lagrangian in terms of the two sets of fermions 4 /j (the 
light quarks) and (the heavy quarks): 


^ = Y^OiO^d, - + i^(^ §z4^/ - S^)), 

f ^ ^ 


where, 




2N 1 ^ _ 1 ^ _ 1 

- b E + y-,Yoi,nf - ^-d,s<^]- 


iV2 _ 


(46) 


(47) 


/=i 


/=! 


We can further refine Eq. (1171) by taking A“ equal to the real part of the right hand side which leads to the additional 
constraint on the field dt /: 


d,{Y^fe^f) = Q. 

/=i 

In Eq. ((151) the term L '2 is the pure gluon Lagrangian independent of the gauge coupling constant: 


(48) 


4 = 


( 49 ) 
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Note that in Eq. (|T7l) the presence of the scalars 5'“ (a = 1...8) is purely optional as there are many ways in which 
one can implement this constraint. Also if are considered large one can integrate out the fermions this leading 
to an equivalent Lagrangian expressed only in terms of the light degrees of freedom f. 


V. DISCUSSION 

When one constructs an effective theory for the low energy degrees of freedom of QCD, be these fermions or hadrons, 
one uses often as a guiding principle the approximate symmetries already established experimentally such as the chiral 
SU{3)l X SU{3)r symmetry or the 17(1) a axial one. In the present work we considered a different approach; thus 
instead of dwelling on symmetries we focused on the intrinsic dynamics encapsulated in the partition function of the 
QCD Lagrangian. Using alternative sets of variables we were able to reproduce the exact partition function obtained 
after integrating out the quark degrees of freedom. It turns out that a partition function depending on a set of fermion 
variables can be described alternatively only by new sets of fermion variables, consider them copies of the first ones, 
that however have a completely different Lagrangian. 

The final result (see Eq. (1461) 1 is an effective Lagrangian with unusual terms depending on two groups of fermions 
(where T/ are the three light quarks and are the three heavy ones) and with couplings which go up to 8 fermion 
interaction terms. This Lagrangian gives an adequate partition function corresponding to that of the original QCD 
Lagrangian for the case when the coupling constant g is large. Moreover it has an intrinsic SU{3)l x SU{3)ii chiral 
symmetry in the two sets of fermion sectors for the case when the quark masses ruf and Mf are set to zero and also 
an SU(S)v symmetry for the case when the two sets of masses are equal within one set. 

The purely quark Lagrangain we obtained can be further processed by using QCD rules to extract the bound states 
of mesons or baryons and their interactions. However the overall method goes far beyond this. One could at any 
intermediate stage introduce directly instead of currents the meson states. In this case one must take into account 
that initially there are 72 degrees of freedom for the light quarks corresponding to a product of 3 colors, 3 flavors and 
8 space time coordinates of an off-shell fermion. This means that the partition function must accommodate only that 
number of mesons, be they scalars, pseudoscalar, vectors etc. whose total number of degrees of freedom sum up 72. 
All alternatives that conveniently express the initial Lagrangian should be considered and all Lagrangians that can 
be constructed in this way are possible outcomes. This is also true for the Lagrangian in Eq. (1461) as it is only one of 
the many possible choices compatible with the initial QCD partition function. 

In this paper we propose a new theoretical effective fermion Lagrangian obtained from QCD not by integrating out 
the gluons but by eliminating these as a result of mocking up the exact partition function of QCD. We shall leave the 
discussion of phenomenological implications of our model for further work. 
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